Optimal metastability-containing sorting via parallel prefix computation

ABSTRACT

In order to provide smaller, faster and less error-prone circuits for sorting possibly metastable inputs, a novel sorting circuit is provided. According to the invention, the circuit is metastability-containing.

RELATED APPLICATIONS

This application is a 371 of International application PCT/EP2019/079861, filed 31 Oct. 2019, the entire contents of which are hereby fully incorporated herein by reference for all purposes. PCT/EP2019/079861 claims the priority benefit of European patent application EP 18000850.0, filed 31 Oct. 2018, the entire contents of which are hereby fully incorporated herein by reference for all purposes.

STATEMENT REGARDING SPONSORED RESEARCH

The project leading to this application has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 716562).

The present invention relates to metastability-containing circuits for sorting an arbitrary number of inputs.

INTRODUCTION

Metastability is a fundamental obstacle when crossing clock domains, potentially resulting in soft errors with critical consequence. As it has been shown that metastability cannot be avoided deterministically, synchronizers are employed to reduce the error probability to tolerable levels. This approach trades precious time for reliability: the more time is allocated for metastability resolution, the smaller the probability of metastability-induced faults.

Recently, a different approach has been proposed, coined metastability-containing (MC) circuits (S. Friedrichs, M. Függer and C. Lenzen, “Metastability-Containing Circuits,” in IEEE Transactions on Computers, vol. 67, no. 8, pp. 1167-1183, 1 Aug. 2018). It accepts a limited amount of metastability in the input to a digital circuit and ensures limited metastability of its output, so that the result is still useful. In particular when sorting inputs, metastability can be contained when sorting inputs arising from time-to-digital converters, i.e., measurement values can be correctly sorted without resolving metastability using synchronizers first.

RELATED WORK

Sorting Networks: Sorting networks sort n inputs from a totally ordered universe by feeding them into n parallel wires that are connected by 2-sort elements, i.e., subcircuits sorting two inputs; these can act in parallel whenever they do not depend on each other's output. A correct sorting network sorts all possible inputs, i.e., the wires are labeled 1 to n such that the i^(th) wire outputs the i^(th) element of the sorted list of inputs. The size of a sorting network is its number of 2-sort elements and its depth is the maximum number of 2-sort elements an input may pass through until reaching the output.

Parallel Prefix Computation: Ladner and Fischer (R. E. Ladner, M. J. Fischer, “Parallel prefix computation”, JACM, vol. 27, no. 4, pp. 831-838, 1980) studied the parallel application of an associative operator to all prefixes of an input string of length l (over an arbitrary alphabet). They give parallel prefix computation (PPC) circuits of depth O(log l) and size O(l) (given a constant-size circuit implementing the operator). A number of additional constructions have been developed for adders, and special cases of the construction by Ladner and Fischer were discovered (in all likelihood) independently, cf. [24]. However, no other construction simultaneously achieves asymptotically optimal depth and size.

It is an object of the invention to provide smaller, faster and less error-prone circuits for sorting possibly metastable inputs.

This object is achieved by a circuit according to independent claim 1. Advantageous embodiments are defined in the dependent claims.

According to an aspect of the invention, CMOS implementations of basic gates realize Kleene logic. The task of comparing inputs can be decomposed into performing a four-valued comparison on each prefix pair of two input strings, followed by inferring the corresponding output bits. Plugging the resulting 2-sort(B) circuits for B-bit inputs into a sorting network for n values readily yields an MC sorting circuit for n valid strings.

The above reduces the task of MC sorting to a parallel prefix computation (PPC) problem, for which circuits that are simultaneously (asymptotically) optimal in depth and size are known due to a celebrated result by Ladner and Fischer (Richard E Ladner and Michael J Fischer. Parallel Prefix Computation. JACM, 27(4):831-838, 1980). According to an aspect of the invention, the inventive circuits can be derived using their framework, which allows for a trade-off between depth and size of the 2-sort circuit. Most prominently, optimizing for depth reduces the depth of the circuit to optimal ┌log B┐, at the expense of increasing the size by a factor of up to 2. However, relying on the construction from Ladner at al. as-is results in a very large fan-out. In a further aspect, the invention proposes reducing fan-out to any number f≥3 without affecting depth, increasing the size by a factor of only 1+O(1/f) (plus at most 3B/2 buffers). In particular, our results imply that the depth of an MC sorting circuit can match the delay of a non-containing circuit, while maintaining constant fan-out and a constant-factor size overhead.

Post-layout area and delay of the designed circuits compare favorably with a baseline provided by a straightforward non-containing implementation.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 shows standard transistor-level implementations of inverter (left), NAND (center), and NOR (right) gates in CMOS technology. The latter can be turned into AND and OR, respectively, by appending an inverter.

FIG. 2 shows finite state machine determining which of two Gray code inputs g, h∈B^(B) is larger. In each step, the machine receives g_(i)h_(i) as input. State encoding is given in square brackets.

FIG. 3 shows An example for a computation of the 2-sort(9) circuit arising from the inventive construction for fan-out f=3. The inputs are g=101010110 and h=101M100000; see Table 10 (FIG. 22 ) for s_(M) ^(i)(g,h) and the output. We labeled each ⋄M by its output. Buffers and duplicated gates (here the one computing 0M) reduce fan-out, but do not affect the computation. Grey boxes indicate recursive steps of the PPC construction; see also FIG. 7 for a larger PPC circuit using the one here in its “right” top-level recursion. For better readability, wires not taking part in a recursive step are dashed or dotted.

FIGS. 4A-4C show the recursion tree T4 (center). Right nodes are depicted black, left nodes gray and leafs are depicted white. The recursive patterns applied at left and right nodes are shown on the left and right, respectively. At the root and its left child, we have that B=B/2; for other nodes, B gets halved for each step further down the tree (where the leaves simply wire their single input to their single output). The left pattern comes in different variants. The basic construction does not incorporate the gray buffers; these will be needed in Section 5.2 to reduce fan-out. The gray wire with index B+1 is present only if B is odd; this never occurs in PPC(C,T_(b)), but becomes relevant when initially applying the left pattern exclusively for k∈N steps (see Theorem 5.8), reducing the size of the resulting circuit at the expense of increasing its depth by k.

FIG. 5 shows comparison of the balanced recursion from [19] and ours. The curves for unbounded fan-out are the exact sizes obtained, whereas “upper bound” refers to the bound from Corollary 5.7; the fan-out 3 curves show that the unbalanced strategy performs better also for the construction from Theorem 5.18 (for f=3 and k=0) we derive next.

FIGS. 6A and 6B show construction of PPC(C,T₄)′. On the left (FIG. 6A), we see the recursion tree, with the aggregation trees separated and shown at the bottom. Inputs are depicted as black triangles. On the right (FIG. 6B), the application of the recursive patterns at the children of the root is shown. Parts 600, 602 will be duplicated in the second step of the construction that achieves constant fan-out; this will also necessitate duplicating some gates in the aggregation trees.

FIG. 7 shows PPC(³⁾(C,T₄). Right recursion steps R_(r) are marked with dark gray, left recursion steps with light gray. The steps at the root (above) and aggregation trees (below) are not marked explicitly. Duplicated gates are depicted in a layered fashion. Dashed lines indicate that a wire is not participating in a recursive step.

FIG. 8 shows a dependence of the size of the modified construction on f. For comparison, the upper bound from Corollary 5.7 on the circuit with unbounded fan-out is shown as well.

FIG. 9 shows an XMUX circuit according to an embodiment of the invention, used to implement ⋄_(M) and out_(M).

FIG. 10 shows constructing 2-sort(B) from out_(M) and PPCM(B−1).

FIG. 11 shows an excerpt from a simulation for 4-bit inputs, where X=M. The rows show (from top to bottom) the inputs g and h, both outputs of the simple non-containing circuit, and both outputs of our design. Inputs g and h are randomly generated valid strings. Columns 1 and 3 show that the simpler design fails to implement a 2-sort(4) circuit.

FIG. 12 shows a comparison of the inventive solution PPC Sort to a standard non-containing one. For the latter, the unexpected delay reduction at B=16 is the result of automatic optimization with more powerful gates, which the inventive solution does not use.

FIGS. 13 to 16 depict Tables 1 to 4, respectively;

FIG. 17 depicts Tables 5a and 5b; and

FIGS. 18 to 22 depict Tables 6 to 10, respectively.

DETAILED DESCRIPTION

We set [N]:={0, . . . , N−1} for N∈

and [i, j]={i, i+1, . . . , j} for i, j∈

, i≤j. We denote

:={0, 1} and

_(M):={0,1, M}. For a B-bit string g∈

_(M) and i∈[1,B], denote by g_(i) its i-th bit, i.e., g=g₁g₂ . . . g_(B). We use the shorthand g_(i,j):=g_(i) . . . g_(j), where i, j∈[1,B] and i≤j. Let par(g) denote the parity of g∈

^(B), i.e, par(g)=Σ_(i=1) ^(B)g_(i) mod 2. For a function f and a set A, we abbreviate f(A):={f(y)|y∈A}.

A standard binary representation of inputs is unsuitable: uncertainty of the input values may be arbitrarily amplified by the encoding. E.g. representing a value unknown to be 11 or 12, which are encoded as 1011 resp. 1100, would result in the bit string 1MMM, i.e., a string that is metastable in every position that differs for both strings. However, 1MMM may represent any number in the interval from 8 to 15, amplifying the initial uncertainty of being in the interval from 11 to 12. An encoding that does not lose precision for consecutive values is Gray code.

A B-bit binary reflected Gray code, rg_(B):[N]→

, is defined recursively. For simplicity (and without loss of generality) we set N:=2^(B). A 1-bit code is given by rg₁(0)=0 and rg₁(1)=1. For B>1, we start with the first bit fixed to 0 and counting with rg_(B−1)(⋅) (for the first 2^(B−1) codewords), then toggle the first bit to 1, and finally “count down” rg_(B−1)(⋅) while fixing the first bit again, cf. Table 1 (FIG. 13 ). Formally, this yields for x∈[N]

$\begin{matrix} {{{{r{g_{B}(x)}}:} = \left\{ \begin{matrix} {0r{g_{B - 1}(x)}} & {{{if}\mspace{14mu} x} \in \left\lbrack 2^{B - 1} \right\rbrack} \\ {1\ r{g_{B - 1}\left( {2^{B} - 1 - x} \right)}} & {{{if}\mspace{14mu} x} \in {\left\lbrack 2^{B} \right\rbrack \smallsetminus \left\lbrack 2^{B - 1} \right\rbrack}} \end{matrix} \right.}.} & (x) \end{matrix}$

As each B-bit string is a codeword, the code is a bijection and the encoding function also defines the decoding function. Denote by

⋅

:

^(B)→[N] the decoding function of a Gray code string, i.e., for x∈[N],

rg_(B)(x)

=x.

For two binary reflected Gray code strings g, h∈

^(B), we define their maximum and minimum as

$\left( {{\max^{rg}\left\{ {g,h} \right\}},{\min^{rg}\left\{ {g,h} \right\}}} \right):=\left\{ \begin{matrix} \left( {g,h} \right) & {{{if}\mspace{14mu}\left\langle g \right\rangle} \geq \left\langle h \right\rangle} \\ \left( {h,g} \right) & {{{if}\mspace{14mu}\left\langle g \right\rangle} < \left\langle h \right\rangle} \end{matrix} \right.$

For example: max^(rg){0011,0100}=max^(rg) {rg _(B)(2),rg _(B)(7)}=0100, min^(rg){0111,0101}=min^(rg) {rg _(B)(9),rg _(B)(10)}=0111.

Inputs to the sorting circuit may have some metastable bits, which means that the respective signals behave out-of-spec from the perspective of Boolean logic. However, they are valid strings in the sense of the invention. Valid strings have at most one metastable bit. If this bit resolves to either 0 or 1, the resulting string encodes either x or x+1 for some x, cf. Table 2 (FIG. 14 ).

More formally, if B∈

and N=2^(B), the set of valid strings of length B is defined as

${S_{rg}^{B}:} = {{r{g_{B}\left( \lbrack N\rbrack \right)}}\bigcup{\bigcup\limits_{x \in {\lbrack{N - 1}\rbrack}}\left\{ {r{g_{B}(x)}*r{g_{B}\left( {x + 1} \right)}} \right\}}}$

The operator * is called the superposition and is defined as

${\forall{i \in {\left\{ {1,\ldots\mspace{14mu},B} \right\}\left( {x*y} \right)_{i}}}}:=\left\{ \begin{matrix} x_{i} & {{{if}\mspace{14mu} x_{i}} = y_{i}} \\ M & {{else}.} \end{matrix} \right.$

The specification of max^(rg) and min^(rg) may be extended to valid strings in the above sense by taking all possible resolutions of metastable bits into account. More particularly, in order to extend the specification of max^(rg) and min^(rg) to valid strings, the metastable closure (Stephan Friedrichs, Matthias Függer, and Christoph Lenzen. Metastability-Containing Circuits. Transactions on Computers, 67, 2018) is used. The metastable closure of an operator on binary inputs extends it to inputs that may contain metastable bits, by considering all possible stable resolutions of the inputs, applying the operator and taking the superposition of the results.

The closure is the best one can achieve w.r.t. containing metastability with clocked logic using standard registers, i.e., when f_(M)(x)_(i)=M, no such implementation can guarantee that the i^(th) output stabilizes in a timely fashion.

If one wants to construct a circuit computing the maximum and minimum of two valid strings, allowing to build sorting networks for valid strings, one also needs to answer the question what it means to ask for the maximum or minimum of valid strings. To this end, suppose a valid string is rg_(B)(x)*rg_(B)(x+1) for some x∈[N−1], i.e., the string contains a metastable bit that makes it uncertain whether the represented value is x or x+1. If one waits for metastability to resolve, the string will stabilize to either rg_(B)(x) or rg_(B)(x+1). Accordingly, it makes sense to consider rg_(B)(x)*rg_(B)(x+1) “in between” rg_(B)(x) and rg_(B)(x+1), resulting in the following total order on valid strings (cf. Table 2, FIG. 14 ).

Definition (<). A total order < is defined on valid strings as follows. For g, h∈

^(B), g<h⇔

g

<

h

. For each x∈[N−1], we define rg_(B)(x)<rg_(B)(x)*rg_(B)(x+1)<rg_(B)(x+1). We extend the resulting relation on S_(rg) ^(B)×S_(rg) ^(B) to a total order by taking the transitive closure. Note that this also defines ≤, via g≤h⇔(g=h∨g<h).

We intend to sort with respect to this order. It turns out that implementing a 2-sort circuit w.r.t. this order amounts to implementing the metastable closure of max^(rg) and min^(rg). In other words, max_(M) ^(rg) and min_(M) ^(rg) are the max and min operators w.r.t. the total order on valid strings shown in Table 2 (FIG. 14 ), e.g., max_(M) ^(rg){1001,1000}=rg ₄(15)=1000, max_(M) ^(rg){0M10,0010}=rg ₄(3)*rg ₄(4)=0M10, and max_(M) ^(rg){0M10,0110}=rg ₄(4)=0110.

Hence, our task is to implement max_(M) ^(rg) and min_(M) ^(rg).

Definition (2-sort(B)). For B∈

, a 2-sort(B) circuit is specified as follows.

Input: g, h∈S_(rg) ^(B)

Output: g′, h′∈S_(rg) ^(B)

Functionality: g′=max_(M) ^(rg){g, h}, h′=min_(M) ^(rg){g, h}.

FIG. 1 shows standard transistor-level implementations of inverter (left), NAND (center), and NOR (right) gates in CMOS technology. The latter can be turned into AND and OR, respectively, by appending an inverter.

The invention seeks to use standard components and combinational logic only. In particular, the behavior of basic gates on metastable inputs may be specified via the metastable closure of their behavior on binary inputs, cf. Table 3 (FIG. 15 ), using the standard notational convention that a+b=OR_(M) (a,b) and ab=AND_(M)(a,b). In this logic, most familiar identities hold: AND and OR are associative, commutative, and distributive, and DeMorgan's laws hold. However, naturally the law of the excluded middle becomes void. For instance, in general, OR(x, x)≠1, as OR(M,M)=M.

It can be shown that the basic CMOS gates shown in FIG. 1 behave according to this logic, i.e. that they implement the truth tables given in Table 3 (FIG. 15 ), thereby justifying the model.

FIG. 2 shows finite state machine determining which of two Gray code inputs g, h∈B^(B) is larger. In each step, the machine receives g_(i)h_(i) as input. State encoding is given in square brackets.

More particularly, FIG. 2 depicts a finite state machine performing a four-valued comparison of two Gray code strings. In each step of processing inputs g, h∈BB, it is fed the pair of ith input bits g_(i)h_(i). In the following, we denote by s^((i))(g,h) the state of the machine after i steps, where s⁽⁰⁾(g,h):=00 is the starting state. For ease of notation, we will omit the arguments g and h of s^((i)) whenever they are clear from context.

Because the parity keeps track of whether the remaining bits are compared w.r.t. the standard or “reflected” order, the state machine performs the comparison correctly w.r.t. the meaning of the states indicted in FIG. 2 .

For all i∈[1,B], we have that max

${{{}_{}^{}\left\{ {g,h} \right\}_{}^{}}{\min^{rg}\left\{ {g,h} \right\}_{i}}} = {{out}\mspace{14mu}\left( {{\overset{l - 1}{\underset{j = 1}{♦}}g_{j}h_{j}},{g_{i}h_{i}}} \right)}$

In order to extend this approach to potentially metastable inputs, all involved operators are replaced by their metastable closure: for i∈[1, B] (i) compute s^((i)), (ii) determine max^(rg){g, h}_(i) and min^(rg){g, h}_(i) according to Table 4 (FIG. 16 ), and finally (iii) exploit associativity of the operator computing the state s(i) in the PPC framework. Thus, we only need to implement ⋄_(M) and the out_(M) (both of constant size), plug them into the framework, and immediately obtain an efficient circuit.

The reader may raise the question why we compute s_(M) ⁽¹⁾ for all i∈[09,B−1] instead of computing only s_(M) ^((B))) with a simple tree of ⋄_(M) elements, which would yield a smaller circuit. Since s_(M) ^((B)) is the result of the comparison of the entire strings, it could be used to compute all outputs, i.e., we could compute the output by out(s_(M) ^((B)), g_(i)h_(i)) instead of out(s^((i−1)) ^(M) , g_(i)h_(i)). However, in case of metastability, this may lead to incorrect results: e.g., for g=0M1 and h=001, we have that s_(M) ⁽³⁾=00*01=0M and out_(M)(0M; g₂h₂)=MM, yet min_(M) ^(rg){g, h}₂=0 (see Tables 6 and 7 (FIGS. 18 and 19 )).

While it is not obvious that this approach yields correct outputs, it may be formally proven that: (P1) ⋄_(M) is associative. (P2) repeated application of ⋄_(M) computes s_(M) ^((i)). (P3) applying out_(M) to s_(M) ^((i−1)) and g_(i)h_(i) results for all valid strings in max_(M) ^(rg){g, h}_(i) min_(M) ^(rg){g, h}_(i). This yields the desired correctness. Regarding the first point, we note the statement that ⋄_(M) is associative does not depend on B. In other words, it can be verified by checking for all possible x, y, z∈B_(;M) ² whether (x⋄_(M)y)⋄_(M)z=x⋄_(M)(y⋄_(M)z).

While it is tractable to manually verify all 3⁶=729 cases (exploiting various symmetries and other properties of the operator), it is tedious and prone to errors. Instead, it was verified that both evaluation orders result in the same outcome by a short computer program, proving the desired associativity of the operator.

For the convenience of the reader, Table 6 (FIG. 18 ) gives the truth table of ⋄_(M). It can be shown that repeated application of this operator to the input pairs g_(j)h_(j)j∈[1, i], actually results in s_(M) ^((i)). This is closely related to the elegant recursive structure of Binary Reflected Gray Code, leading to the important observation that if in a valid string there is a metastable bit at position m, then the remaining B−m following bits are the maximum code word of a (B−m)-bit code.

It may be observed that for g∈S_(rg) ^(B), if there is an index 1≤m<B such that g_(m)=M then g_(m+1,B)=10^(B−m−1).

The reasoning is based on distinguishing two main cases: one is that s_(M) ^((i)) contains at most one metastable bit, the other that s_(M) ^((i′))MM. Each of these cases can be proven by technical statements.

It may further be observed that if |res(s_(M) ^((i)))|≤2 for any i∈[B+1], then res(s_(m) ^((i)))=⋄_(j=1) ^(i)res(g_(j)h_(j)).

The operator out: B²×B²→B² is the operator given in Table 4 (FIG. 16 ) computing max^(rg){g, h}_(i) min^(rg){g, h}_(i) out of s^((i−1)) and g_(i)h_(i). For convenience of the reader, we provide the truth table of out_(M) in Table 7 (FIG. 19 ).

FIG. 3 shows an example for a computation of the 2-sort(9) circuit arising from the inventive construction for fan-out f=3. The inputs are g=101010110 and h=101M100000; see Table 10 (FIG. 22 ) for s_(M) ^(i)(g,h) and the output. More particularly, table 10 (FIG. 22 ) shows an example run of the FSM in FIG. 2 on inputs g=101010110 and h=101M100000. We drop s_(M) ⁽⁹⁾s(9)M, as it is not needed to compute g′₉h′₉. We labeled each ⋄M by its output. Buffers and duplicated gates (here the one computing 0M) reduce fan-out, but do not affect the computation. Grey boxes indicate recursive steps of the PPC construction; see also FIG. 7 for a larger PPC circuit using the one here in its “right” top-level recursion. For better readability, wires not taking part in a recursive step are dashed or dotted.

In order to derive a small circuit from the above, we make use of the PPC framework by Ladner and Fischer. They described a generic method that is applicable to any finite state machine translating a sequence of B input symbols to B output symbols, to obtain circuits of size O(B) and depth O(log B). They reduce the problem to a parallel prefix computation (PPC) task by observing that each input symbol defines a restricted transition function, whose compositions evaluated on the starting state yield the state of the machine after the corresponding number of steps. This matches our needs, as we need to determine s_(M) ^((i)) for each i∈[B]. However, their generic construction involves large constants. Fortunately, we have established that ⋄_(M):B_(M) ²×B_(M) ² is an associative operator, permitting us to directly apply the circuit templates for associative operators they provide for computing s_(M) ^((i))=(⋄_(M))_(j=1) ^(i)g_(j)h_(j) for all i∈[B]. Accordingly, only these templates are discussed.

We revisit the part of the framework relevant to our construction, also providing a minor improvement on their results in the process. To this end, we first formally specify the PPC task for the special case of associative operators.

Definition 5.1 (PPC_(⊕)(B)). For associative ⊕: D×D→D and B∈N, a PPC_(⊕)(B) circuit is specified as follows.

Input: d∈D^(B),

Output: π∈D^(B),

Functionality: π_(i)=⊕_(j=1) ^(i)d_(j) for all i∈[1, B].

In our case, ⊕=⋄_(M) and D=B_(M) ² the method by Ladner et al. provides a family of recursive constructions of PPC_(⊕) circuits. They are obtained by combining two different recursive patterns.

FIG. 4 shows the recursion tree T4 (center). Right nodes are depicted black, left nodes gray and leaves are depicted white. The recursive patterns applied at left and right nodes are shown on the left and right, respectively. At the root and its left child, we have that B=B/2; for other nodes, B gets halved for each step further down the tree (where the leaves simply wire their single input to their single output). The left pattern comes in different variants. The basic construction does not incorporate the gray buffers; these will be needed to reduce fan-out. The gray wire with index B+1 is present only if B is odd; this never occurs in PPC(C,T_(b)), but becomes relevant when initially applying the left pattern exclusively for k∈N steps, reducing the size of the resulting circuit at the expense of increasing its depth by k.

More particularly, suppose that C and P are circuits implementing ⊕ and PPC_(⊕)(┌B/2┐) for some B∈N, respectively. Then applying the recursive pattern given at the left of FIG. 4 (i) with B:=B and without the rightmost gray line if B is even and (ii) with B:=B−1 if B is odd yields a PPC_(⊕)(B) circuit. It has depth 2d(C)+d(P) and size at most (B−1)|C|+|P|. Moreover, the last output is at depth at most d(C)+d(P) of the circuit.

The second recursive pattern, shown in FIG. 4 c , avoids to increase the depth of the circuit beyond the necessary d(C) for each level of recursion. Assume for now that B is a power of 2. We represent the recursion as a tree T_(b), where b:=log B, given in the center of FIG. 6A. It has depth b with all leafs in this depth, and there are two types of nonleaf nodes: right nodes (filled in black) have two children, a left and a right node, whereas left nodes (filled in gray) have a single child, which is a right node. T_(b) is essentially a Fibonacci tree in disguise.

Definition. T₀ is a single leaf. T₁ consists of the (right) root and two attached leaves. For b≥2, T_(b) can be constructed from T_(b−1) and T_(b−2) by taking a (right) root r, attaching the root of T_(b−1) as its right child, a new left node l as the left child of r, and then attaching the root of T_(b−2) as (only) child of l.

The recursive construction is now defined as follows. A right node applies the pattern given in FIG. 4 to the right, where R_(l) is the circuit (recursively) defined by the subtree rooted at the left child, R_(r) is the circuit (recursively) defined by the subtree rooted at the right child, and B=2^(b−d−1), where d∈[b] is the depth of the node. A left child applies the pattern on the left, where the recursively used circuit R_(c) is defined by the subtree rooted at its child and B=2^(b−d), where d∈[b] is the depth of the node. The base case for a single input and output is simply a wire connecting the input to the output, for both patterns. As b=log B and each recursive step cuts the number of inputs and outputs in half, the base case applies if and only if the node is a leaf. Note that the figure shows the recursive patterns at the root and its left child, where B=2^(b−1) is always even (i.e., in this recursive pattern, the gray wire with index B+1 is never present); when applying the patterns to nodes further down the tree, B and B are scaled down by a factor of 2 for every step towards the leaves.

In the following, denote by PPC(C,T_(b)) the circuit that results from applying the recursive construction described above to the base circuit C implementing ⊕. Moreover, we refer to the i^(th) input and output of the subcircuit corresponding to node ν∈T_(b) as d_(i) ^(ν) and π_(i) ^(ν), respectively.

It may be shown that If C implements ⊕, PPC(C,T_(b)) is a PPC_(⊕)(2^(b))circuit and PPC(C,T_(b)) has depth b·d(C).

It remains to bound the size of the circuit. Denote by F_(i), i∈N, the i^(th) Fibonacci number, i.e., F₁=F₂=1 and F_(i+1)=F_(i)+F_(i−1) for all 2≤i∈N. Then it may be shown that PPC(C,T_(b)) has size (2^(b+2)−F_(b+5)+1)|C|.

Asymptotically, the subtractive term of F_(b+5) is negligible, as F_(b+5)∈(1/√{square root over (5)}+0(1))((1+√{square root over (5)})/2)^(b+5)⊆O(1.62^(b)); however, unless B is large, the difference is substantial. We also get a simple upper bound for arbitrary values of B. To this end, we “split” in the recursion such that the left branch is “complete,” while applying the same splitting strategy on the right. This is where our construction differs from and improves on the method of Ladner et al. They perform a balanced split and obtain an upper bound of 4B on the circuit size.

It follows that for B∈N and circuit C implementing ⊕, set b:=┌log B┐. Then a PPC_(⊕)(B) of depth ┌log B┐d(C) and size smaller than (5B−2^(b)−F_(b+3))|C|≤(4B−F_(b+3)) exists.

We remark that one can give more precise bounds by making case distinctions regarding the right recursion, which for the sake of brevity we omit here. Instead, we computed the exact numbers for B≤70.

FIG. 5 shows comparison of the balanced recursion from Ladner et al and ours. The curves for unbounded fan-out are the exact sizes obtained, whereas “upper bound” refers to the above-given bound; the fan-out 3 curves show that the unbalanced strategy performs better also for the construction (for f=3 and k=0) we derive next.

The construction derived from iterative application of the above results can be combined with PPC(C,T_(b)), achieving the following trade-off; note that if B=2^(b) for b∈N, then F_(┌log B┐−k+3) can be replaced by F_(b−k+5).

Suppose C implements ⊕. For all k∈[┌log B┐+1] and b∈N, there is a PPC_(⊕)(B) circuit of depth (┌log B┐+k)d(C) and size at most

$\left( {{\left( {2 + \frac{1}{2^{k - 1}}} \right)B} - F_{{\lceil{\log\; B}\rceil} - k + 3}} \right){C}$

FIGS. 6A and 6B show the construction of PPC(C,T₄)′. On the left (FIG. 6A), we see the recursion tree, with the aggregation trees separated and shown at the bottom. Inputs are depicted as black triangles. On the right (FIG. 6B), the application of the recursive patterns at the children of the root is shown. Parts 600, 602 will be duplicated in the second step of the construction that achieves constant fan-out; this will also necessitate duplicating some gates in the aggregation trees.

The optimal depth construction incurs an excessively large fan-out of Θ(B), as the last output of left recursive calls needs to drive all the copies of C that combine it with each of the corresponding right call's outputs. This entails that, despite its lower depth, it will not result in circuits of smaller physical delay than simply recursively applying the construction from FIG. 4A. Naturally, one can insert buffer trees to ensure a constant fan-out (and thus constantly bounded ratio between delay and depth), but this increases the depth to Θ(log² B+d(C)log B).

We now modify the recursive construction to ensure a constant fan-out, at the expense of a limited increase in size of the circuit. The result is the first construction that has size O(B), optimal depth, and constant fan-out.

In the following, we denote by f≥3 the maximum fanout we are trying to achieve, where we assume that gates or memory cells providing the input to the circuit do not need to drive any other components. For simplicity, we consider C to be a single gate.

We proceed in two steps. First, we insert 2B buffers into the circuit, ensuring that the fan-out is bounded by 2 everywhere except at the gate providing the last output of each subcircuit corresponding to a right node. In the second step, we will resolve this by duplicating such gates sufficiently often, recursively propagating the changes down the tree. Neither of these changes will affect the output of the circuit or its depth, so the main challenges are to show our claim on the fan-out and bounding the size of the final circuit.

Step 1: Almost Bounding Fan-Out by 2

Before proceeding to the construction in detail, we need some structural insight on the circuit.

For node ν∈T_(b), define its range R_(υ) and left-count α_(υ) recursively as follows.

-   -   If υ is the root, then R_(υ)=[1,2^(b)] and α_(υ)=0.     -   If υ is the left child of p with R_(p)=[i, i+j], then R_(υ)=[i,         i+(j+1)/2] and α_(υ)=α_(p).     -   If υ is the right child of right node p with R_(p)=[i, i+j],         then R_(υ)=[i+(j+1)/2+1, i+j] and α_(υ)=α_(p).     -   If υ is the right child of left node p, then R_(υ)=R_(p) and         α_(υ)=α_(p)+1.

Suppose the subcircuit of PPC(C,T_(b)) represented by node ν∈T_(b) in depth d∈[b+1] has range R_(υ)=[i, i+j].

-   -   Then     -   (i) it has 2^(b−d) inputs,     -   (ii) j=2^(b−d+α) ^(ν) −1,     -   (iii) if υ is a right node, all its inputs are outputs of its         childrens' subcircuits, and     -   (iv) if υ is a left node or leaf, only its even inputs are         provided by its child (if it has one) and for odd k∈[1,2^(b−d)],         we have that

d_(k)^(v) = ⊕_(k^(′) = i + (k − 1))^(i + k2^(α_(v) − 1))2^(α_(v))d_(k^(′)).

This leads to an alternative representation of the circuit PPC(C,T_(b)), see FIGS. 6A-6B, in which we separate gates in the recursive pattern from FIG. 4A that occur before the subcircuit R_(c). Adding the buffers we need in our construction, this results in the modified patterns given in FIG. 6B. The separated gates appear at the bottom of FIG. 6A: for each leaf υ of T_(b), there is a tree of depth α_(υ) aggregating all of the circuit's inputs from its range. Each non-root node in an aggregation tree provides its output to its parent. In addition, one of the two children of an inner node in the tree must provide its output as an input to one of the subcircuits corresponding to a node of T_(b′), cf. Property (iv) above.

From this representation, we will derive that the following modifications of PPC(C,T_(b)) result in a PPC_(⊕)(2^(b)) circuit PPC(C,T_(b))′, for which a fan-out larger than 2 exclusively occurs on the last outputs of subcircuits corresponding to nodes of T_(b).

-   -   1) Add a buffer on each wire connecting a non-root node of any         of the aggregation trees to its corresponding subcircuit (see         FIG. 6A).     -   2) For the subcircuit corresponding to left node l with range         R_(l)=[i, i+j], add for each even k≤j (i.e., each even k but the         maximum of j+1) a buffer before output π_(k) ^(l) (see bottom of         FIG. 6B).     -   3) For each right node r with range [i, i+j], add a buffer         before output π_((j+1)/2) ^(r) (see top of FIG. 6B).

With the exception of gates providing the last output of subcircuits corresponding to nodes of T_(b) (600, 602 in FIG. 6B), fan-out of PPC(C,T_(b))′ is 2. Buffers or gates driving an output of the circuit drive nothing else.

It remains to count the inserted buffers. The following helper statement will be useful for this, but also later on.

Denote by L_(b)⊆T_(b) the set of leaves of T_(b). Then |L_(b)|=F_(b+2) and Σ_(ν∈L) _(b) 2^(α) ^(ν) =2^(b).

Next, consider the recurrence given by L′₀=1, L′₁=2, and L′_(b)=L′_(b−1)+2L′_(b−2) for b≥2; the factor of 2 assigns twice the weight to the subtree rooted at the child of the root's left child, thereby ensuring that each leaf is accounted for with weight 2^(α) ^(ν) . This recurrence has solution 2^(b).

Denote by s the size of a buffer. Then |PPC(C,T _(b))′|=|PPC(C,T _(b))|+(2^(b)+2^(b−1) −F _(b+3))s.

Step 2: Bounding Fan-Out by f

In the second step, we need to resolve the issue of high fan-out of the last output of each recursively used sub circuit in PPC(C,T_(b))′. Our approach is straightforward. Starting at the root of T_(b) and progressing downwards, we label each node υ with a value α(υ) that specifies a sufficient number of additional copies of the last output of the sub circuit represented by υ to avoid fan-out larger than f. At right nodes, this is achieved by duplicating the gate computing this output sufficiently often, 600, 602 in FIG. 6B (top). For left nodes, we simply require the same number of duplicates to be provided by the sub circuit represented by their child (i.e., we duplicate the 602 wire in the bottom recursive pattern shown in FIG. 6B). Finally, for leaves, we will require a sufficient number of duplicates of the root of their aggregation tree; this, in turn, may require to make duplicates of their descendants in the aggregation tree.

We start by defining α(υ) and then argue how to use these values for modifying the circuit to obtain our fan-out f circuit. Afterwards, we will analyze the increase in size of the circuit compared to PPC(C,T_(b))′.

Definition (α(υ)). Fix b∈N₀. For ν∈T_(b) in depth d∈[b+1], define

${{a(v)}:} = \left\{ \begin{matrix} 0 & {{if}\mspace{14mu} v\mspace{14mu}{is}\mspace{14mu}{the}\mspace{14mu}{root}} \\ \frac{{a(p)} + 2^{b - d}}{f} & {{if}\mspace{14mu} v\mspace{14mu}{is}\mspace{14mu}{the}\mspace{14mu}{left}\mspace{14mu}{child}\mspace{14mu}{of}\mspace{14mu} p} \\ \frac{a(p)}{f} & {{if}\mspace{14mu} v\mspace{14mu}{is}\mspace{14mu}{the}\mspace{14mu}{right}\mspace{14mu}{child}\mspace{14mu}{of}\mspace{14mu}{right}\mspace{14mu}{node}\mspace{14mu} p} \\ p & {{if}\mspace{14mu} v\mspace{14mu}{is}\mspace{14mu}{the}\mspace{14mu}({only})\mspace{14mu}{child}\mspace{14mu}{of}\mspace{14mu}{left}\mspace{14mu}{node}\mspace{14mu}{p.}} \end{matrix} \right.$

Suppose that for each leaf ν∈T_(b), there are └a(ν)┘ additional copies of the root of the aggregation tree, and for each right node ν∈T_(b), we add └a(ν)┘ gates that compute (copies of) the last output of their corresponding sub circuit of PPC(C,T_(b))′. Then we can wire the circuit such that all gates that are not in aggregation trees have fan-out at most f, and each output of the circuit is driven by a gate or buffer driving only this output.

It remains to modify the aggregation trees so that sufficiently many copies of the roots' output values are available.

Consider an aggregation tree corresponding to leaf ν∈T_(b) and fix f≥3. We can modify it such that the fan-out of all its non-root nodes becomes at most f, there are └a(ν)┘ additional gates computing the same output as the root, and at most (fa(ν))/(f−2)+(2^(α) ^(ν) ⁻¹)/(f−1) gates are added.

Finally, we need to count the total number of gates we add when implementing these modifications to the circuit.

For f≥3, define PPC^((f))(C,T_(b)) by modifying PPC(C,T_(b))′ according to Lemmas 5.15 and 5.16. Then, with λ₁:=(1+√{square root over (5)})/4, |PPC^((f))(C,T_(b))| is bounded by

${{{PPC}\left( {C,T_{b}} \right)}^{\prime}} + {2^{b}\left( {\frac{1}{{2f} - 2} + \frac{2}{f - 2} + {O\left( \frac{\lambda_{1}^{b}}{f^{2}} \right)}} \right){{C}.}}$

We summarize our findings in the following:

Suppose that C implements ⊕, buffers have size s and depth at most d(C), and set λ₁:=(1+√{square root over (5)})/4. Then for all k∈[b+1], b∈N₀ and f≥3, there is a PPC_(⊕)(2^(b)) circuit of fan-out f, depth (b+k)d(C), and size at most

${\left( {2^{b + 1} + {2^{b - k}\left( {2 + \frac{{5f} - 6}{{2f^{2}} - {6f} + 4} + {O\left( \frac{\lambda_{1}^{b}}{f^{2}} \right)}} \right)}} \right){C}} + {\left( {2^{b} + 2^{b - k - 1}} \right){s.}}$

Due to space constraints, we refrain from analyzing the size of the construction for values of B that are not powers of 2. However, in FIG. 8 we plot the exact bounds (without buffers) for k=0 and selected values of f against B.

FIG. 7 shows, as an example for the overall resulting construction, PPC⁽³⁾(C,T₄). Right recursion steps R_(r) are marked with dark gray, left recursion steps with light gray. The steps at the root (above) and aggregation trees (below) are not marked explicitly. Duplicated gates are depicted in a layered fashion. Dashed lines indicate that a wire is not participating in a recursive step.

FIG. 8 shows a dependence of the size of the modified construction on f. For comparison, the upper bound on the circuit with unbounded fan-out is shown as well.

FIG. 9 shows an XMUX circuit according to an embodiment of the invention, used to implement ⋄_(M) and out_(M).

First, we need to specify implementations of the sub circuits computing ⋄_(M) and out_(M).

From Tables 5a and 5b (FIG. 17 ), for s, b∈B² we can extract the Boolean formulas (s⋄b)₁ =s ₁ s ₂ +s ₁ b ₁ +s ₂ b ₁ (s⋄b)₂ =s ₁ s ₂ +s ₁ b ₂ +s ₂ b ₂ out(s,b)₁ =s ₁ b ₂ +s ₂ b ₁ +b ₁ b ₂ out(s,b)₂ =s ₁ b ₂ +s ₂ b ₁ +b ₁ b2.

In general, realizing a Boolean formula f by replacing negation, multiplication, and addition by inverters, AND, and OR gates, respectively, does not result in a circuit implementing f_(M) ¹ However, we can easily verify that the above formulas are disjunctions of all prime implicants of their respective functions. As one can manually verify, these formulas evaluate to the truth tables given in Tables 6 and 7 (FIGS. 18 and 19 ) and in this special case the resulting circuits do implement the closure—provided the gates behave as in Table 3 (FIG. 15 ), which the implementations given in FIG. 2 do. Using distributive laws (recall that these also hold in Kleene logic), the above formulas can be rewritten as (s⋄b)₁ =s ₁( s ₂ +b ₁)+ s ₂ b ₁ (s⋄b)₂ =s ₂( s ₁ +b ₂)+ s ₁ b ₂ out(s,b)₁ =b ₁(b ₂ +s ₂)+b ₂ s ₁ out(s,b)₂ =b ₂(b ₁ +s ₁)+b ₁ s ₂. ¹ For instance, (s⋄b)₁=s₁b₁+s ₂b₁ as Boolean formula, but the two expressions differ when evaluated on s₁=s ₂=1 and b₁=M. The circuits resulting from the different formulas are implementations of a multiplexer (with select bit b1) and its closure, respectively.

We see that, in fact, a single circuit with suitably wired (and possibly negated) inputs can implement all four operations. As for sel₁=sel₂ the circuit implements a multiplexer with select bit sel₁, we refer to it as extended multiplexer, or XMUX for short. Its functionality is specified by XMUX(sel₁,sel₂ ,x,y):=y(x+sel₂)+x sel₁.

Table 8 (FIG. 20 ) lists how to map inputs to compute ⋄_(M) and out_(M).

We note that this circuit is not a particularly efficient XMUX implementation; a transistor-level implementation would be much smaller. However, our goal here is to verify correctness and give some initial indication of the size of the resulting circuits—a fully optimized ASIC circuit is beyond the scope of this article. The size of the implementation may be slightly reduced by moving negations. Due to space limitations, we refrain from detailing this modification here, but note that FIG. 12 and table 9 (FIG. 21 ) consider it.

We now have all the pieces in place to assemble a containing 2-sort(B) circuit. As stated above, ⋄_(M) is associative. Thus, from a given implementation of ⋄_(M) (e.g., two copies of the circuit from FIG. 9 with appropriate wiring and negation, cf. Table 8 (FIG. 20 )) we can construct PPC_(⋄) _(M) (B−1) circuits of small depth and size, as shown above. We can combine such a circuit with an out_(M) implementation (again, two XMUXes with appropriate wiring and negation will do) to obtain our 2-sort(B) circuit.

FIG. 10 shows constructing 2-sort(B) from out_(M) and PPCM(B−1).

The correctness of this construction follows from the above explanations, where we can plug in any PPC_(⋄) _(M) (B−1) circuit. For the circuits derived by relying on the XMUX circuit from FIG. 9 , we independently confirmed this via simulation.

More particularly, we implemented the design given in FIG. 10 on register transfer-level using the PPC_(⋄) _(M) (B−1) circuit described above for k=0.³ Quartus by Altera was used for design entry, which in our case mainly consists of checking correct implementation. After design entry, we used ModelSim by Altera for behavioral simulation. Note that we must not simulate the preprocessed Quartus output, because processing may compromise metastability-containing behavior. Instead, we simulate pure VHDL. Metastable signals are simulated using VHDL signal X, because its behavior matches the worst-case behavior assumed for M.

FIG. 11 shows an excerpt from a simulation for 4-bit inputs, where X=M. The rows show (from top to bottom) the inputs g and h, both outputs of the simple non-containing circuit, and both outputs of our design. Inputs g and h are randomly generated valid strings. Columns 1 and 3 show that the simpler design fails to implement a 2-sort(4) circuit.

For the implementation of PPC_(⋄) _(M) (B−1) we used the basic circuits, i.e., we did not make use of the extension to constant fan-out. We exhaustively checked the design from FIG. 10 for B up to 12 (and all k accordingly). Simulation shows that the design works correct for several levels of recursion, e.g., when regarding B=1 and B=2 as simple base cases, B=12 implies 3 levels of recursion for both patterns. We refrained from simulating the constant fan-out construction, because it simply repeats replicates intermediate results without adding functionality.

After behavioral simulation we continue with a comparison of our design and a standard sorting approach Bin-comp(B). As mentioned earlier, the 2-sort(B) implementation given in FIG. 9 is slightly optimized by pulling out a negation from the operators in every recursive step [3]. After design entry as described above we use Encounter RTL Compiler for synthesis and Encounter for place and route. Both tools are part of the Cadence tool set and in both steps we use NanGate 45 nm Open Cell Library as a standard cell library.

Since metastability-containing circuits may include additional gates that are not required in traditional Boolean logic, Boolean optimization may compromise metastability containing properties. Accordingly, we were forced to disable optimization during synthesis of the circuits.

FIG. 12 shows a comparison of the inventive solution PPC Sort to a standard non-containing one. For the latter, the unexpected delay reduction at B=16 is the result of automatic optimization with more powerful gates, which the inventive solution does not use.

As a binary benchmark Bin-comp was used: In short, Bin-comp consists of a simple VHDL statement comparing two binary encoded inputs and outputting the maximum and the minimum, accordingly. It follows the same design process as 2-sort, but then undergoes optimization using a more powerful set of basic gates. For example, the standard cell library provides prebuilt multiplexers. These multiplexers are used by Bin-comp, but not by 2-sort. We stress that these more powerful gates provide optimized implementations of multiple Boolean functions, yet each of them is still counted as a single gate. Thus, comparing our design to the binary design in terms of gate count, area, and delay disfavors our solution. Moreover, we noticed that the optimization routine switches to employing more powerful gates when going from B=8 to B=16 (cf. FIG. 12 ), resulting in a decrease of the delay of the Bin-comp implementation.

Nonetheless, our design performs comparably to the non-containing binary design in terms of delay, cf. FIG. 12 and Table 9 (FIG. 21 ). This is quite notable, as further optimization of our design is possible by optimizing it on the transistor level, with significant expected gains. The same applies to gate count and area, where a notable gap remains. Recall, however, that the Bin-comp design hides complexity by using more advanced gates and does not contain metastability.

We emphasize that we refrained from optimizing the design by making use of all available gates or devising transistor-level implementations, since such an approach is tied to the utilized library or requires design of standard cells.

In conclusion, we demonstrated that efficient metastability containing sorting circuits are possible. Our results indicate that optimized implementations can achieve the same delay as non-containing solutions, without a dramatic increase in circuit size. This is of high interest to an intended application motivating us to design MC sorting circuits: fault tolerant high-frequency clock synchronization. Sorting is a key step in envisioned implementations of the Lynch-Welch algorithm with improved precision of synchronization. The complete elimination of synchronizer delay is possible due to the efficient MC sorting networks presented in this article; enabling an increment of the rate at which clock corrections are applied, significantly reducing the negative impact of phase drift of local clock sources on the precision of the algorithm.

More generally speaking, MC circuits like those presented here are of interest in mixed signal control loops whose performance depends on very short response times. When analog control is not desirable, traditional solutions incur synchronizer delay before being able to react to any input change. Using MC logic saves the time for synchronization, while metastability of the output corresponds to the initial uncertainty of the measurement; thus, the same quality of the computational result can be achieved in shorter time. Note that our circuits are purely combinational, so they can be used in both clocked and asynchronous control logic.

Examples of such control loops are clock synchronization circuits, but MC has been shown to be useful for adaptive voltage control and fast routing with an acceptable low probability of data corruption as well. 

The invention claimed is:
 1. A metastability-containing sorting circuit, the sorting circuit comprising a sorting network for sorting a number n of input strings, each string having length B, wherein the input strings are Gray coded, wherein the sorting network comprises at least two 2-sort circuits C for sorting each pair of input strings, and wherein a size of the sorting circuit is within the order of the size of the input strings O(B) and a maximum fanout f≥3 of the sorting circuit is constant, characterized in that a depth of the sorting circuit is within the order of log(B).
 2. The circuit of claim 1, wherein the Gray code is a binary reflected Gray code.
 3. The circuit of claim 1, wherein buffers are used to bound the fanout.
 4. A transistor-level implementation of the circuit of claim
 1. 5. The circuit of claim 1, wherein a depth of the sorting circuit is (┌log(B)┐+k)d(C), where k∈{0, . . . , ┌log (B)┐} and d(C) is a depth of a 2-sort circuit. 